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Abstract 



Orthogonality criterion is used to shown in a very simple and general way that 
anomalous bound-state solutions for the Coulomb potential (hydrino states) do 
not exist as bona fide solutions of the Schrodinger, Klein-Gordon and Dirac equa- 
tions. 



An alleged tightly-bound state of hydrogen with strong singularity of the 
eigenfunction at the origin (called a hydrino state) has received considerable at- 
tention in the literature [Ij. The order of magnitude of the atomic size (Bohr 
radius) as well as the energy of the hydrogen atom in its ground state just derived 
from the Heisenberg uncertainty principle, even in a relativistic framework, should 
be enough to disqualify hydrino states. However, in a recent Letter, Dombey pQ 
rejects the solution of the three-dimensional Klein-Gordon equation, previously 
derived by Naudts [2], as well as the solution of the two-dimensional Dirac equa- 
tion, by resorting to a few fair arguments. In addition, Dombey presents a solid 
argument founded on the Hermiticity of the Hamiltonian for the Klein-Gordon 
case and a suggestion of similar treatment for the three-dimensional Dirac case. 
In the wake of Dombey's suggestion, this Letter presents such a general criterion 
for banishing hydrino states in the context of the standard quantum mechanics. 

The time-independent Schrodinger equation 
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and the time-independent Klein-Gordon equation 

- h 2 c 2 V 2 ^ + M 2 cV = (E — Vf i> (2) 
with spherically symmetric potentials admit eigenfunctions in the form 

i , = ^l Ylm {9^) (3) 
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where k denotes the principal quantum number plus other possible quantum 
numbers, u k is a square-integrable function (J °° dr \u k \ 2 = 1) and Y™ are the 
orthonormalized spherical harmonics (J dflY^Yj^ = S^Smfh), with I = 0, 1, 2, . . . 
and m = —I, —/ + 1, in such a way that 

H eti u k = (E ett ) k u k (4) 
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. &_JL + v ^— + n2 ' (/+1) for the Klein-Gordon equation 
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and 

{E for the Schrodinger equation 

(6) 
E 2Mc 2<: ^ or ^ e Klein-Gordon equation 
Meanwhile, the time- independent Dirac equation is 

Hij k = E k ij k} H = a-p + (3Mc 2 + V (7) 

where «j and j3 in the standard (or Dirac-Pauli) representation are given by the 
2x2 block matrix form 

«=[ ° "A P=(l ° 1, i = l,2,3 



CTj J \ -1 

and <Tj are the 2x2 Pauli matrices. Its eigenfunction has a spinorial structure 
given by [3] 



where k = =F (j + 1/2), with the minus sign for aligned spin (j = 1 + 1/2) 
and the plus sign for unaligned spin (j = 1 — 1/2). Here, y* mj (with j = 
1/2, 3/2, 5/2, . . . and rrij = — j, — j + 1, . . . , j) are the orthonormalized spinor 
spherical harmonics resulting from the coupling of two-dimensional spinors to the 
eigenstates of orbital angular momentum. The normalization of the Dirac spinor 
(Jo°° ^ r (l/fc| 2 + \dk\ 2 ) — 1) requires that f k and g k are square-integrable functions. 
Using the identity a ■ V = a ■ r — CT and the properties a • fy* = — 3^-^ 
and ex • Z} 1 ! = — (« + 1) X K ™ i there results that one can write 

J ''*7 v f J llt j 

H eii $k = E k $ k (10) 

where 

fhc(-f + Z) V-Mc 2 \ , ( f k \ 



It is instructive to examine the radial solutions in the neighbourhood of the 
origin for the Coulomb potential —hca/r (a is the coupling constant) because 
u, f and g must behave better than r -1 / 2 at the origin in order to guarantee 
their square integrability. An estimative of the asymptotic behaviour of the ra- 
dial solutions for small r can be obtained by neglecting the terms of order r~ n 
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compared with the terms of order r~( n+l \ where n = 0, 1, 2, in the Schrodinger, 
Klein-Gordon and Dirac differential equations. 

As r — > the terms behaving as r~ 2 (r _1 ) dominate in the Schrodinger case 
for / ^ (/ = 0) in such a way that 

d 2 v 

r 2 ^-l(l + l)u = 0, for/ t^O (12) 
dr l 

2 d 2 u 2Mca n . , n 

r 2 — H — ru = 0, for Z = (13) 
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Substituting the Frobenius series expansion 

oo 

u = J2W +n , a ^0 (14) 

n=0 

into (IT21 and ffl~3]) yields the quadratic indicial equation (obtained when n = 0): 
z/ (z/ — 1) = I (/ + 1), which has the solutions v\ = I + 1 and z/ 2 = — I. It happens 
that there is no recurrence formula for I ^ (a n = for n ^ 0) so that the 
general solution for ( 1121) can be written as u = v4;r i+1 + Bir~ l . On the other 
hand, for I = the power series expansion fails with v — v-i because only the 
first coefficient of the series can be defined. It means that the Frobenius method 
furnishes only one solution for I = 0. A second solution linearly independent 
of Ui = a^r + the solution related to v — Vi, can be found by writing u 2 = 

oo 

aui In (r) + b n r U2+n with a ^ and bo ^ 0. It follows that the general solution 

n=0 

for I = can be written as u = A (r - ^fr 2 H )+B (l - ^ff^r In (r) + ■ ■ •) . 

The above results imply that only choosing Bi = B = gives a behaviour 
at the origin which can lead to square-integrable solutions. For short, all the 
normalizable radial solutions (u/r) are regular at the origin and behave as r l as 
r — ► 0, for all I. 

As for the Klein-Gordon case, the terms behaving as r -2 , for all values of I, 
outweighs the other terms in such a manner that one can write 

r 2 ^- [l(l + l)-a 2 ]u = (15) 

for small r. The power series expansion ( fl4|) gives the indicial equation v [y — 1) = 
I (I + 1) — a 2 which has the roots 

u± =}-±J(l + ^) -a 2 , a<l + \ (16) 
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The argumentation given below is based on (I16p and can be better understood 
by observing Fig. [TJ where v + and z/_, as a function of a, are plotted on the 
same grid. The constraint v > —1/2 which ensures the square integrability of u 
demands that for a < 1/2 there are just a S-wave (I = 0) solution for v = V- 
although all the values of I are allowed for v = v + . Since v < 1 for a < 1/2, 
both S-wave solutions, in the sense of u/r, diverge at the origin. As a increases 
some normalizable solutions diverging at the origin for v = v_ become possible 
whereas some normalizable solutions for v = v + become divergent. In any case, 
as the coupling constant takes critical values (a c = l c + 1/2, l c — 0, 1, 2, . . .) only 
the solutions for / > l c are allowed. 

In analogy with the solutions in the neighbourhood of the origin of the Schrodinger 
and Klein-Gordon equations, we look for solutions of the asymptotic Dirac equa- 
tion 



df , . 
r— + k} -ag 
dr 



dg , 
r—-Kg + af 
dr 
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in the form of power series 



n=0 



n=0 



Now the coefficient of r v (n = 0) gives the system of indicial equations 



[y + h) a — abo 



aa + [v — k) bo 







(19) 



and the two possible values of v: 

v± = 



a' 



a < \k\ 



(20) 



Fig. [2] illustrates (120]) as a function of the coupling constant for the lowest values 
of |«|. The requirement of normalizability [y > —1/2) implies that for a < y/3/2 
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only the positive root is allowed and the radial solutions (f/r and g/r) are found 
for all the values of \k\ (\k\ = 1,2,3,...), although the radial solutions with \k\ = 1 
diverge at the origin. The negative root allows divergent radial solutions with 
\k\ — 1 for y/3/2 < a < 1. Nevertheless, there are no normalizable solutions with 
\k\ — 1 for a > 1. Continuing this process, as the coupling constant increases 
starting with a = 1, some normalizable solutions for v = v + turn out to be 
divergent whereas some normalizable solutions diverging at the origin for v — V- 
become possible. As a takes critical values (a = \k\ c , \k\ c = 1, 2, 3, . . .) only the 
solutions for \k\ > \k\ c are allowed. This means that the radial solutions for the 
smallest values of \k\ become disallowed one after the other as a increases. 

In the standard quantum mechanics an observable such as the energy is repre- 
sented by an Hermitian operator, whose set of eigenf unctions constitutes a basis 
so that every arbitrary wave function can be expanded in one and only one way 
in terms of the eigenfunctions. Besides square-integrability, appropriate bound- 
ary conditions must be imposed on the eigenfunctions of an eigenvalue problem. 
Square-integrability requires that the eigenfunction vanishes at the infinity and 
the boundary condition at the origin for a singular potential, as the Coulomb 
potential, comes naturally into existence by demanding that the Hamiltonian is 
Hermitian, viz. 



where ipf. is an eigenfunction corresponding to an eigenvalue Ek- In passing, note 
that a necessary consequence of Eq. ( |2~TT) is that the eigenfunctions corresponding 
to distinct effective eigenvalues are orthogonal. Identifying H e±f with H and 
E eii with E in (12TI) . integrating by parts and recalling the orthonormality of 
the spherical harmonics, it is easy to show that Uf-(r) for the Schrodinger and 
Klein-Gordon cases must satisfy the following constraint [l]-[5] 



For the Dirac case, identifying H ef± with H in (121]) . one finds a constraint involv- 
ing the upper and lower components, namely 



For the Schrodinger case, the square-integrable solution for u (u— > r l+1 as 
r — » 0) satisfies (|22|) and its corresponding radial solution u/r is regular at the 




(21) 




(22) 




(23) 
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origin. 

For the Klein-Gordon case, only the solutions for u less singular than r 1//2 
can satisfy the orthogonality criterion. One sees, then, that the square-integrable 
solution must behave as r u+ , as r — > 0, where v + is one of the two possibilities 
of dH]), corresponding to the thick lines in Fig. HJ Therefore, square-integrable 
solutions of the Klein-Gordon equation satisfying the orthogonality criterion are 
such that for a < 1/2 all the values of I are allowed. The I = radial solution, 
in the sense of u/r, has a divergence at the origin. As a increases some regular 
normalizable solutions become divergent. In any case, as the coupling constant 
takes critical values (a c — l c + 1/2, l c = 0, 1, 2, . . .) we are left with the solutions 
for I > l c . 

For the Dirac case, only the regular solutions for / and g can satisfy the 
orthogonality criterion. This means that the square-integrable solutions must 
behave at the origin as r u+ , where u + is the positive square root of fT20l) . corre- 
sponding to the thick lines in Fig. [2j Here, square-integrable solutions satisfying 
the orthogonality criterion are such that, for a < 1, the radial solutions (f/r and 
g/r) are found for all the values of \k\ (\k\ = 1, 2, 3, . . .), with the radial solutions 
with \k\ — 1 diverging at the origin. Nevertheless, there are no normalizable 
solutions with \k\ = 1 for a > 1. Continuing this process, as a increases starting 
with a = 1, some solutions turn out to be divergent. As a takes critical values 
(a = \k\ c , \k\ c = 1,2,3,...) only the solutions for \k\ > \k\ c are allowed. This 
means that the radial solutions for the smallest values of \k\ become disallowed 
one after the other as a increases. 

In summary, a very simple and general criterion has been presented to reject 
hydrino states in the context of the standard quantum mechanics. Square inte- 
grability is sufficient enough to exclude singular wave functions in the Schrodinger 
equation, but not in the Klein-Gordon and Dirac ones. For the relativistic equa- 
tions with the Coulomb potential, singular wave functions are allowed. The or- 
thogonality criterion, though, imposes an additional constraint in such way that 
the would-be relativistic square-integrable solutions for hydrino states, related to 
the thin lines in Fig. [1] and Fig. [21 are not acceptable. Therefore, only radial 
solutions behaving at the origin as r _1 / 2+e , with e > 0, for the Klein-Gordon case, 
and r _1+e for the Dirac case, are physically acceptable solutions. 



6 



Acknowledgments 

The author is indebted to an anonymous referee for very constructive remarks 
and to CNPq and FAPESP for partial financial support. 



References 

[1] N. Dombey, Phys. Lett. A 360 (2006) 62 and references therein. 
[2] J. Naudts, |physics/05071~93j 

[3] W. Greiner, Relativistic Quantum Mechanics, Wave Equations, Springer- 
Verlag, Berlin, 1990. 

[4] K. M. Case, Phys. Rev. 80 (1950) 797. 

[5] D. Xianxi et al., Phys. Rev. A 55 (1997) 2617 and references therein. 



7 




Figure 1: The roots of the indicial equation, v + (thick line) and z/_ (thin line), 
as a function of the coupling constant for the Klein-Gordon case (a for I = 0, 
b for I = 1, c for I = 2 and d for / = 3). The shaded area represents the zone 
corresponding to radial solutions diverging at the origin. 
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Figure 2: The solutions of the system of indicial equations, u + (thick line) and 
v_ (thin line), as a function of the coupling constant for the Dirac case (a for 
\k\ = 1, b for \k\ = 2 and c for |ac| = 3). The shaded area represents the zone 
corresponding to radial solutions diverging at the origin. 
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